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1.  Introduction 


SASPRO  II,  an  acronym  for  5pare  and  Server  Provisioning  (model 
two),  is  a versatile  FORTRAN  package  that  gives  provisioning  levels  of 
spares  inventory  and  repair  capacity  required  to  support  a population  of 
randomly  failing  items  which,  upon  failure,  are  (1)  dispatched  to  the  re- 
pair facility,  and  (2)  replaced  by  a spare  if  one  is  available.  This  pa- 
per describes  in  detail  the  problem  environment,  the  program  options,  the 
input  required  to  run  the  program,  and  the  output  provided  by  the  pro- 
gram. Sample  runs  are  also  shown  for  each  of  the  program  options. 

2.  Problem  Environment 


Consider  a population  of  items  containing  certain  significant 
parts;  for  example,  a fleet  of  aircraft  containing  key  avionics  gear,  a 
fleet  of  ships  with  modular  engine  components,  or  a group  of  milling  ma- 
chines, where  the  entire  machine  itself  is  the  key  "part."  These  "parts" 
randomly  fail  and  require  repair.  Spare  parts  are  also  needed  so  that  upon 
failure,  the  spare  can  be  utilized  to  replace  the  failed  part  and  the  item 
put  back  into  service.  It  is  desired  to  determine  how  many  spares  and  how 
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many  repaii  channels  are  required  to  support  the  system  at  a desired  ser- 
vice level  while  minimizing  costs. 

The  system  is  shown  schematically  in  Figure  1.  We  consider  only 
a single  part- type  at  a time,  which  has  its  own  spares  pool  and  dedicated 
repair  channels.  For  example,  for  a fleet  of  gas  turbine  propelled  ships, 
the  gas  turbine  engine  has  two  components — a gas  generator  and  a power 
turbine.  Each  must  have  dedicated  repair  channels  and  its  own  spares 
pool.  Thus  SASPRO  II  would  treat  each  component  in  turn,  being  utilized 
to  provision  first  for  a population  of  gas  generators  and  then  for  a pop- 
ulation of  power  turbines. 

When  a unit  in  the  operating  population  fails,  a spare  is  requested 
at  the  same  time  the  unit  is  dispatched  for  repair.  If  a spare  is  not 
available,  the  request  is  backlogged  and  units  coming  out  of  repair  are 
used  in  removing  the  backlog.  When  there  is  no  backlog  of  requests  for 
spares,  units  coming  out  of  repair  go  into  the  spares  inventory.  Repair 
times  as  well  as  failure  times  are  treated  as  random  variables  and  with 
the  proper  assumptions  (to  be  mentioned  below),  this  stochastic  process 
can  be  readily  modeled  as  a finite  source  queueing  system,  often  referred 
to  as  "the  machine  repair  problem  with  spares,"  which,  in  addition,  also 
fits  a two-stage  cyclic  queueing  model.  Thus  SASPRO  IT.  uses  a standard 
queueing  model  for  the  stochastic  process  [see  GROSS,  KAHN,  and  MARSH 
(1977)]. 

The  assumptions  required  for  using  SASPRO  II  are  that  times  to 
failure  and  repair  times  are  exponentially  distributed  random  variables. 
These  assumptions  allow  the  employment  of  the  standard  finite  source 
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queueing  theory  to  determine  probabilities  of  various  numbers  of  units 
in  repair  at  any  given  time.  From  this,  system  service  levels  (number  of 
units  operating,  availability  of  spares,  etc.)  can  readily  be  computed. 

In  order  to  achieve  a specified  service  level,  certain  combinations  of 
spares  and  servers  (repair  channels)  are  required.  Using  costs  associ- 
ated with  purchasing  and  holding  spares,  and  costs  associated  with  build- 
ing and  operating  repair  channels,  SASPRO  II,  through  a heuristic  optimi- 
zation routine,  finds  the  "best"  combination  of  spares  and  servers  that 
meets  the  service  level  constraint.  While  the  optimization  algorithm  is 
heuristic,  it  does  guarantee  a feasible  solution  and  also  calculates  the 
exact  cost  of  the  solution  so  that  the  user  can  "manually"  perturb  the 
heuristic  solution  to  find  better  ones  if  they  exist. 

3.  Modes  of  Operation 

The  program  has  two  modes  of  operation,  dynamic  and  static.  The 
former  is  advised  for  initial  year  provisioning  when  population  sizes, 
failure  rates,  and  repair  times  may  be  changing  significantly.  Popula- 
tion size  changes  may  occur  because  units  are  put  into  operation  gradu- 
ally, thus  building  up  to  a full  strength  population  over  a period  of 
several  years.  For  example,  it  was  anticipated  to  build  a fleet  of  256 
gas  turbine  powered  ships,  starting  in  the  first  year  with  ten  ships  and 
building  up  to  full  strength  over  a ten-year  period.  Because  of  new 
technology,  engines  on  ships  introduced  in  the  later  years  were  expected 
to  have  smaller  failure  rates  (be  more  reliable)  while  due  to  learning, 
repair  times  were  also  expected  to  be  smaller  in  the  later  years. 

In  designing  support  systems  for  which  it  is  necessary  to  deter- 
mine the  number  and  location  of  depots  the  static  mode  is  useful  [see 
GROSS  and  PINKUS  (1978)].  In  this  situation  the  population  is  at  full 
strength,  technological  advances  and  learning  are  complete,  and  condi- 
tions are  very  close  to  static.  Running  tiroes  and  input  requirements 
are  greatly  teduced  when  operating  in  the  static  node. 

The  dynamic  mode  allows  for  changing  population  sizes,  failure 
rates,  and  repair  times,  as  well  as  for  changing  costs  on  a year  by  year 
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basis.  A set  of  input  information  is  required  for  each  year  in  the  plan- 
ning horizon.  An  item  repaired  and  placed  back  into  spares  inventory  (or 
operation)  in  year  i is  assumed  to  have  the  same  failure  rate  as  a new 
item  introduced  in  year  i. 

All  costs  are  turned  into  an  equivalent  beginning  of  year  payment 
and  then  discounted  according  to  where  the  year  is  in  the  planning  hori- 
zon, so  that  at  any  year  i the  program  gives  the  present  worth  of  the 
cumulative  sum  of  the  discounted  costs  up  to  and  including  year  i.  The 
final  value  for  the  last  year  is  then  the  present  worth  of  the  sum  of  dis- 
counted costs  over  the  entire  planning  horizon. 

The  costs  that  are  considered  in  SASPRO  II  are  purchase  costs; 
salvage  values;  and  annual  holding  and  operating  costs  associated  with 
spares  and  repair  channels,  respectively;  unit  transportation  and  repair 
costs;  and  component  improvement  program  (CIP)  investment  costs.  Purchase 
costs  and  salvage  values  are  in  dollars  per  spare  or  repair  channel.  Op- 
erating costs  of  each  channel  and  holding  costs  for  spares  are  in  dollars 
per  year  per  spare  or  repair  channel.  Transportation  and  repair  costs  are 
in  dollars  per  unit  per  repaired  item  and  component  improvement  program 
costs  are  in  dollars  per  year. 

In  determining  the  discounted  algorithm  costs,  salvage  values  are 
not  realizable  until  the  end  of  the  planning  horizon,  even  though  spares 
or  repair  channels  are  retired  prior  to  that;  also,  operating  costs  and 
holding  costs  are  assumed  to  be  incurred  every  year  until  the  end  of  the 
planning  horizon,  even  if  spares  or  channels  are  retired  earlier.  Trans- 
portation repair  and  CIP  costs  are  not  explicitly  considered  by  the  algo- 
rithm. 


The  assumptions  that  salvage  values  are  not  received  until  the  end 
of  the  planning  horizon  and  that  operating  or  holding  costs  are  not  re- 
duced when  spares  or  repair  channels  are  retired  early  are  necessitated  by 
the  heuristic  optimization  algorithm  employed.  Further,  if  a spare  or  re- 
pair channel  is  retired  during  the  planning  horizon  and  is  required  again  a 
few  years  later,  it  must  be  repurchased.  Again,  this  assumption  is  re- 
quired because  of  the  nature  of  the  heuristic  optimization  algorithm. 
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Once  the  heuristic  optimal  values  of  spares  and  repair  channels 
are  obtained,  the  true  cost  and  true  present  worth  can  then  be  deter- 
mined. The  true  costs  are  the  actual  annual  costs.  In  calculating 
these,  the  assumptions  on  salvage  values,  and  operating  and  holding  costs 
are  no  longer  required;  that  is,  the  salvage  values  are  received  when- 
ever spares  or  repair  channels  are  retired,  and  the  operating  or  holding 
costs  are  being  incurred  only  for  those  spares  or  repair  channels  actu- 
ally present  each  year.  Transportation,  repair,  and  CIP  costs  are  also 
included.  Ti.  cue  present  worth  at  any  year  i is  the  cumulative  sum 
of  the  discounted  true  costs  up  to  and  including  year  i.  The  final  value 
for  the  last  year  is  then  the  true  present  worth  of  the  sum  of  discounted 
true  costs  over  the  entire  planning  horizon. 

In  the  static  mode  fewer  problems  arise  since  spares  and  channels 
are  not  added  or  retired.  Costs  for  this  mode  of  operation  are  converted 
to  expected  equivalent  end-of-year  payments  over  the  planned  life  of  the 
system. 


Details  of  the  algorithm  cost  calculations  and  the  true  cost  cal- 
culations are  provided  in  Section  10. 

4.  Service  Level  Constraint  Options 

There  are  two  options  available  for  specifying  service  perform- 
ance. The  first,  referred  to  as  spares  availability , sets  a limit  on  the 
percentage  of  requests  for  spares  that  are  met  from  on-shelf  spares  in- 
ventory (also  called  fill  rate);  that  is, 

Number  of  Requests  for  Spares  per  Year  Honored  Immediately  > . . 

Number  of  Requests  for  Spares  per  Year  — ’ ' 1 

where  A is  specified  by  the  user  and  0 < A < 1 . 

The  second  criterion  for  service  performance,  called  "fleet" 
availability , sets  a level  for  the  percentage  of  time  a certain  portion 
of  the  population  desired  to  be  in  operation  is  actually  operating;  that 
is, 

Pr{>  BM  units  are  uo}  > A . (2) 
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Both  3 and  A are  specified  by  the  user  where  0 < 3 < 1 , and  M is 
the  population  (fleet)  size  excluding  spares. 

Suppose,  for  example,  we  wish  to  have  100  machines  in  operation 
(M=100).  Wnen  a machine  fails,  a spare  machine  is  "plugged  in"  if  one  is 
available.  We  might  specify  that  a service  level  constraint  be  (1)  the 
percentage  of  requests  for  spare  machines  filled  immediately  from  on-hand 
spares  is  at  least  90%  (option  1:  A=.9),  or  (2)  at  least  95%  of  the  ma- 
chines are  operating  85%  of  the  time  (option  2:  3=-95,  A=.85). 


5.  Population  Average  Failure  Rate  Options 

Using  the  average  failure  rate  for  each  year  allows  for  the  incor- 
poration of  changing  component  reliability  as  the  years  progress.  There 
are  two  options  available  for  specifying  population  average  failure  rate: 
averaging  the  failure  rates  and  averaging  the  mean  time  between  failures 
(or  removals,  denoted  by  MTBR) . 


Population  average  failure  rate  computed  by  averaging  the  failure 
rates  is  calculated  by  the  formulae  (we  refer  to  this  as  rate  averaging) : 


{<VMi-i>Ai + Ri-xxi-i + <"i-rRi-i)Xi-i) ' Mi  • Mi > Mi-i  <3) 

i=2  3 

< Vixi-1  + <Mi-rVi)xi-i> ' M1-1  ' Mi  £ M1-1 


For  averaging  the  MTBRs,  the  population  average  failure  rate  is 
given  by  (we  refer  to  this  as  time  averaging): 


X1  = A1 


A"1  = I 


{(MrMi-i)/xi  + Vi/xi-i  + (Mi-rRi-i)/Ai-i} ' Mi  * Wi 


(4) 


^Ri_l/Xi_l  + (Mi-l-Ri-l)/Ai-i^  / Mi-l 


1=2,3 


where 
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Mj,  = component  population  size,  year  i 
A^  = component  failure  rate,  year  i 

R,  = expected  number  of  components  repaired,  year  i. 

Which  is  the  better  averaging  method  to  use  depends  on  whether 
there  are  many  or  few  machines  operating  simultaneously.  For  the  many 
machine  case,  rate  averaging  yields  more  accurate  results  while  for  the 
few  machine  case,  time  averaging  is  more  appropriate  [see  GROSS  and  INCE 
(1978)].  Note  that  both  Equations  (3)  and  (A)  assume  that  the  likelihood 
of  a given  item  failing  more  than  once  in  a time  period  is  negligible. 

6.  Perturbation  Options 

The  heuristic  algorithm  operates  on  a year-by-year  basis  and  hence 
has  the  limitation  that  it  does  not  "look  ahead."  Further,  it  treats  op- 
erating costs  and  salvage  values  in  a very  approximate  way  and  does  not 
explicitly  consider  repair,  transportation,  or  CIP  costs  at  all.  But  it 
does  yield  a feasible  solution;  that  is,  one  that  will  meet  the  service 
level  constraint.  When  considering  the  entire  planning  horizon,  the  heu- 
ristic algoiithm  may  (and  probably  does)  not  prove  optimal.  Therefore, 
after  obtaining  the  heuristic  optimal  solution  it  may  be  advisable  to  make 
some  perturbations  by  visually  selecting  other  values.  By  exercising  the 
perturbation  option  the  program  will  use  the  perturbed  solution  values  as 
if  they  were  the  optimal  solution  (without  going  through  the  heuristic  al- 
gorithm) and  will  print  out  the  availabilities  and  true  costs,  allowing 
comparison  to  those  given  by  the  heuristic  optimal  solution. 

7.  Input  Data 

Table  I shows  the  data  that  are  required  as  input  for  SASPRO  II. 

Most  input  parameters  are  self-explanatory  but  a few  require  fur- 
ther comment.  The  A shown  in  Equations  (1)  and  (2)  is  AVL,  while  BETA 
is  the  3 shown  in  Equation  (2). 

When  using  the  algorithm  (KTC  = 1),  initial  values  CO  and  YO  must 
be  read  in  for  number  of  servers  and  spares,  respectively.  In  the  dynam- 
ic mode,  after  the  first  year  the  program  uses  the  previous  year's  values 
for  CO  and  YO  as  initial  values,  thus  the  CO  and  YO  fields  can  be  left 


- 7 - 


T-391 


TABLE  I 

INPUT  REQUIREMENTS 


Variable 

Name 

Description 

Symbol 

on 

Printout 

AVI, 

Desired  Availability 

AVL 

BETA 

Desired  Percent  of  Population  Up 

BETA 

C 

Initial  Value — Number  of  Repair  Channels 

CO 

CIC 

Carrying  Cost  per  Spare  ($/yr/spare) 

CIC 

Cl  PC 

Component  Improvement  Cost  ($/yr) 

CIPC 

CPSER 

Repair  Channel  (Server)  Purchase  Cost  ($/channel) 

CPSER 

CPSP 

Spare  Purchase  Cost  ($/spare) 

CPSP 

H 

Operating  Hours  per  Year  per  Item  (hrs) 

H 

JD 

„ T (=1:  Rate  Average 

Averaging  Rate  Option  Indicator  j=2;  Tlme  Average 

KEYWD 

Mode  Option  Indicator  | 

(=0:  Static 

KTC 

„ . _ (=1:  Heuristic  Opt.  Algorithm 

Perturbation  Indicator  { n . 

(=2:  Perturbation  Option 

KWRITE 

t , , _ - , T , • =0:  Not  Print 

Intermediate  Output  Option  Indicator  , . 

=1:  Print 

KZ 

o . ^ ^ -r  , . (=0:  Spare  Avail 

Service  Criterion  Option  Indicator  { . 

r (=1:  Fleet  Avail 

NY  EARS 

Planning  Horizon  Length  (vrs) 

YRS 

OCPSER 

Operating  Cost  of  a Channel  ($/vr/channel) 

OCPSER 

R 

Yearly  Interest  Rate 

RATE 

RM 

Population  Size 

M 

RMTBR 

Mean  Time  Between  Removals  (hrs) 

MTBR 

ST 

Average  Turn  Around  Time  (days) 

1/MU 

SVPSER 

Salvage  Value  of  a Channel  ($/channel) 

SVPSER 

SVPSP 

Salvage  Value  of  a Spare  ($/spare) 

SVPSP 

URC 

Unit  Repair  Cost  ($/unit) 

URC 

UTC 

Unit  Transportation  Cost  ($/unit) 

UTC 

Y 

Initial  Value — Number  of  Spares 

YO 
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blank  on  the  card  sets  for  every  year  in  the  planning  horizon  after  the 
first.  The  closer  the  initial  values  are  to  the  final  values  (deter- 
mined by  SASPRO  II),  the  fewer  iterations  of  the  heuristic  optimization 
algorithm  are  required.  However,  one  may  use  CO  = YO  = 1 if  so  desired. 

A set  of  cost  inputs  (CIC,  CIPC,  CPSER,  CPSP,  OCPSER,  SVPSER, 

SVPSP,  URC,  UTC)  is  required  for  each  year  in  the  horizon  in  the  dynamic 
mode.  This  allows  one  to  account  for  inflation  and  technological  innova- 
tions. The  Component  Improvement  Program  Cost  (CIPC)  is  the  annual  ex- 
penditure required  to  achieve  a given  MTBR  schedule  (the  MTBR  which  must 
be  inputted  for  each  year  in  the  horizon)  for  the  dynamic  mode  of  opera- 
tion, or  to  maintain  the  MTBR  achieved  when  operating  in  the  static  mode. 

The  MTBR  value  is  the  actual  mean  time  to  failure  of  each  unit 
when  operating  continuously.  If  items  do  not  operate  continuously  but 
are  required  for,  say,  only  H hours  per  year  on  the  average,  the  mean 
failure  rate  actually  used  in  the  queueing  model  portion  of  SASPRO  II  is 
lowered  accordingly.  If  items  do  operate  around  the  clock,  H = 365  * 

24  = 8760.  If,  for  example,  each  unit  in  a population  of  items  has  an 
MTBR  of  1000  hours  but  is  called  upon  to  operate,  on  the  average,  only 
half  the  time  (H  = 4380  hours),  the  effective  MTBR  used  in  the  program 
is  raised  to  2000  hours  (failure  rate  cut  in  half) . The  user  specifies 
H and  MTBR  and  SASPRO  II  automatically  makes  the  adjustment.  The  reader 
is  referred  to  BARZILY,  GROSS,  and  KAHN  (1977)  for  a discussion  of  the 
adequacy  of  this  procedure  to  account  for  noncontinuous  operation.  The 
above  reference  also  discusses  the  SASPRO  II  assumptions,  when  operating 
in  the  dynamic  mode,  that  (1)  the  population  attains  instantaneous  steady- 
state  each  year  at  average  values,  and  (2)  the  population  consists  of  non- 
identical units  (with  respect  to  mean  time  to  failure),  which  are  treated 
as  identical  by  weighted  averaging  according  either  to  Equation  (3)  or 
(4).  Gross  and  Ince  (1978)  further  discuss  this  latter  problem. 

The  parameters  KEYWD,  KZ,  JD,  KTC,  and  KWRITE  are  the  option 
flags.  Setting  KTC  = 2 puts  SASPRO  II  in  the  perturbation  mode;  setting 
KTC  = 1 causes  SASPRO  II  to  operate  with  the  heuristic  optimization  algo- 
rithm. Putting  KEYWD  = 1 sets  SASPRO  II  in  the  dynamic  mode;  setting 
KEYWD  = 0 allows  SASPRO  II  to  operate  in  the  static  mode.  Designating 
KZ  = 1 puts  che  service  level  constraint  on  fleet  availability;  KZ  = 0 
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sets  the  service  level  constraint  on  spares  availability.  A B must  be 
specified  when  KZ  = 1.  Putting  JD  ■ 2 sets  the  population  average  fail- 
ure rate  calculation  to  averaging  MTBRs;  putting  JD  =•  1 averages  failure 
rates.  For  KWRITE  = 0,  intermediate  output  will  also  be  printed;  for 
KWRITE  = 1 only  final  output  is  printed. 

Table  II  and  Figure  2 give  the  card  layout  required  for  the  input 
information.  There  are  eleven  cards  needed  for  the  static  mode  and  nine 
plus  two  cards  for  each  year  in  the  planning  horizon  required  for  dynamic 
mode  operation.  The  input  requirements  for  static  mode  operation  are 
similar  to  those  required  for  a one-year  planning  horizon  dynamic  run. 
However,  the  output  cost  values  given  in  the  static  mode  are  the  expected 
end  of  year  payments  adjusted  over  an  NYEARS  life,  while  the  costs  of  a 
single  year  dynamic  mode  run  are  the  present  worth  of  expenditures  for 
that  year. 

8.  Output  from  SASPRO  II 

SASPRO  II  gives  the  heuristic  optimum  combinat.on  of  spares  and. 
repair  channels  needed  to  meet  the  service  level  constraint  (or  the  ac- 
tual service  level  for  an  inputted  set  of  spares  and  repair  channels)  and 
also  provides  the  costs  associated  with  this  solution.  For  the  static  op- 
eration mode  there  is  a single  line  of  output  with  all  cost  values  being 
the  expected  equivalent  annual  expenditure  over  the  NYEAR  system  life.  For 
the  dynamic  mode  of  operation  there  is  a line  of  output  for  each  year,  the 
costs  outputted  being  the  expected  present  worth  of  the  cumulative  sum  of 
discounted  costs  up  to  and  including  that  year  as  well  as  the  costs  for 
that  particular  year,  as  given  by  both  the  algorithm  and  exact  calculation. 
Also  given  ai  output  are  the  heuristic  optimum  combinations  of  servers  and 
spares  (when  operating  in  the  optimization  mode);  the  average  system  fail- 
ure rate,  which  in  the  static  mode  is  the  same  as  the  inputted  failure  rate 
calculated  from  the  MTBR  and  H values  read  in,  and  in  the  dynamic  mode  is  a 
weighted  average  [according  either  to  Equation  (3)  or  (4)]  of  the  units  in 
the  population  which  were  introduced  and  repaired  in  various  years  at  dif- 
ferent values;  the  average  number  of  units  repaired;  and  the  actual  avail- 
ability achieved  (always  > AVL  when  using  the  heuristic  algorithm).  Another 
output  quantity  shown  is  ASTAR,  the  percentage  of  time  the  population  is 
called  upon  to  operate  (ASTAR  » H/8760);  this  serves  as  a check  on  the  H 
value  put  in.  The  output  quantities  with  definitions  are  shown  in  Table  III. 

- 10  - 


- 


T-391 


TABLE 


II 


CARD  LAYOUT  FOR  INPUT 


Card 

Number 

Input  Data  Parameter(s) 

Format 

Columns 

1 

Title  (any  desired  by  user) 

— 

1-80 

2 

NYEARS 

12 

1-2 

3 

R 

F8.5 

1-8 

4a 

AVL,  BETA 

F8.5,  F8.5 

1-8,  9-16 

5l 

KZ 

12 

1-2 

6C 

KEYWD 

12 

1-2 

7C 

JD 

12 

1-2 

8£ 

KTC 

12 

1-2 

9f 

KWRITE 

12 

1-2 

10 

See  Figure  2:  One  set  required  for 

each  year  in  dynamic  planning  hori- 

“I 

zcn;  one  set  only  for  static  mode. 

For  Spares  Avail  Option,  BETA  may  be  set  at  any  value. 


b _ 0 -*■  Spares  Availability 

1 -*■  Fleet  Availability 


CKEYWD  = Static  Mwd! 

( 1 ■>  Dynamic  Mode 

djD  f 1 -*•  Averaging  Failure  Rates 
j 2 -*•  Averaging  MTBRs 

eKTC  1 Heuristic  Algorithm 
( 2 ->  Perturbation  Option 


fKWRITE 


0 -*■  Not  Print  Intermediate  Output 

1 ■+  Print  Intermediate  Output 
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MTBR  1/KU  H CIPC  CPSER  ' CPSP  I URC 


F5.0  F7.0  F5.0  F5.0 


TABLE  III 


OUTPUT  QUANTITIES 


Name 

Description 

YR 

Actual  year  represented 

M 

Population  size  year  i (from  input) 

FR 

Failure  rate  of  a unit  purchased  or  repaired  in  year  i 
[failures/day  = (1/MTBR)  • (H/8760)  • 24] 

FRBAR 

Average  failure  rate  of  a typical  unit  (failures/day  = 
weighted  average  of  various  units  purchased  or  repaired 
in  all  years  ur  to  and  including  i) 

AS  TAR 

Average  percent  of  time  population  is  called  upon  to  op- 
erate (H/8760) 

C 

Heuristic  optimum  number  of  repair  charnels  required  in 
year  i 

Y 

Heuristic  optimum  number  of  spares  required  in  year  i 

AVAIL 

Availability  achieved 

RBAR 

Average  number  of  units  repaired  in  year  i 

COST 

Costs,  as  considered  by  the  heuristic  algorithm,  expend- 
ed in  year  i dynamic  mode  or  equivalent  yearly  average 
expenditure  in  static  mode 

PR-WORTH 

Present  worth  of  sum  of  discounted  algorithm  costs  up  to 
and  including  year  i,  dynamic  mode;  same  as  COST  for 
static  mode 

TRUE-COST 

True  costs  expended  in  year  i dynamic  mode  or  true 
equivalent  yearly  average  expenditure  in  static  mode 

TRUE-PW 

True  present  worth  of  sum  of  discounted  costs  up  to  and 
including  year  1,  dynamic  mode;  same  as  TRUE-COST  for 
static  mode 
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9.  Sample  Runs 

We  illustrate  the  model  options  by  presenting  eight  sample  runs 
as  given  below  in  Table  IV. 

TABLE  IV 
SAMPLE  HUNS 


Sample 
Run  No 

Options 

Planning 

Horizon 

Mode 

Failure  Rate 
Computation 

Service  Level 
Constraint 

1 

Dynamic 

Heur.  Opt. 

Rate  Avg. 

Fleet  Avail. 

2 

Dynamic 

Heur.  Opt. 

Rate  Avg. 

Spares  Avail. 

3 

Dynamic 

Heur.  Opt. 

Time  Avg. 

Fleet  Avail. 

4 

Dynamic 

Heur.  Opt. 

Time  Avg. 

Spares  Avail. 

5 

Static 

Heur.  Opt. 

Fleet  Avail. 

6 

Static 

Heur.  Opt. 

— 

Spares  Avail. 

7 

Dynamic 

Heur.  Opt. 

Rate  Avg. 

Fleet  Avail. 

8 

Dynamic 

Perturb. 

Rate  Avg. 

Fleet  Avail. 

A listing  of  the  input  cards  for  these  runs  is  given  in  Figure  3. 

For  the  eight  cases,  there  is  a total  of  196  data  input  cards  (6 [9  + 

(2  x 10)1  + 2 [ 11 ] ) . 

The  associated  output  for  the  first  six  cases  (Sample  Runs  1 tc  6) 
is  given  in  Figure  4.  In  Figure  5,  output  for  the  last  two  cases  (Sample 
Runs  7 and  8)  is  shown.  Run  7 is  similar  to  Run  1 but  the  spare  purchase 
cost  in  the  initial  year  is  reduced  from  617  to  400.  The  heuristic  algo- 
rithm solution  does  not  change.  However,  by  inspection  it  seems  that  per- 
turbing the  1979  and  1980  Y values  from  3 and  4,  respectively,  to  6 
should  give  a better  solution  since  the  initial  spare  purchase  cost  is 
relatively  cheap.  Exercising  the  perturbation  mode  (KTC**2)  in  Sample  Run 
8 shows  this  to  be  true  by  comparing  the  TRUE-PW  values  for  the  final  year. 
The  input  requirements  for  Sample  Run  8 would  be  identical  to  those  for 
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Figure  4. — continued 


Figure  4. — continued 


SAMPLE  RUN  7 

OPTIONS  I DYNAMIC,  HEURISTIC  OPT,  RATE  AVC,  FLEET  AVAIL 
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Figure  5. — Illustration  of  perturbation  option. 
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Sample  Run  7 except  KTC  is  set  to  2 and  all  CO,YO  values  must  be  punched 
in  and  would  be  equal  to  the  C,Y  values  of  Sample  Run  7,  except  for 
changing  the  1979  and  1980  YOs  to  6. 

Not*  that  there  is  a large  difference  between  the  algorithm  costs 
and  the  true  costs.  Much  of  this  is  due  to  (URC  + UTC)R  , and  while  R 
is  affected  by  the  choices  of  the  decision  variables  C and  Y,  the  effect 
on  optimality  is  of  a secondary  nature  since  R changes  only  very  slightly 
for  vast  differences  in  C and  Y when  all  else  is  constant  (compare,  for 
example.  Sample  Runs  5 and  6 or  years  1979  and  1 980  in  Sample  Runs  7 and 
8).  While  C1PC  could  have  a sizable  effect  on  optimal  costs  because  it 
is  directly  related  to  attainable  failure  rates,  this  can  be  studied  via 
a sensitivity  type  of  analysis;  that  is,  rerunning  with  a variety  of  CIPC 
programs  and  their  associated  failure  rate  schedules.  In  the  Sample  Runs 
1 through  8,  CIPC  (as  well  as  UTC)  was  set  to  zero.  In  Sample  Run  9, 
shown  in  Figure  6,  we  do  a five-year  dynamic  horizon  with  conditions  the 
same  as  for  Run  1,  except  that  CIPC  is  set  at  400,  5C0,  600,  200,  200  and 
UTC  is  set  at  5,  5,  10,  10,  10,  respectively  over  the  five-year  period. 
Although  the  algorithm  solution  came  out  the  same  as  lor  Run  1,  the  algo- 
rithm costs  differ  somewhat  from  Run  1 due  to  a five-year  rather  than  a 
ten-year  anticipated  horizon.  The  true  cost  and  true  present  worth  differ 
to  account  for  the  added  CIPC  and  UTC  costs. 

10 . Intermediate  Output,  Cost  Functions,  and 

the  Heuristic  Optimization  Algorithm 

Also  provided  (if  so  desired  by  setting  KWRITE-1)  as  output  are 
intermediate  values  of  Y and  C which  "step  up"  from  Y0  and  CO,  showing 
the  operation  of  the  algorithm  at  each  iteration.  Briefly,  the  heuristic 
algorithm  works  as  follows. 

For  the  dynamic  mode,  the  true  present  worth  of  the  sum  of  dis- 
counted yearly  costs  over  a dynamic  horizon  of  K years  is  given  by 
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TRUE-PW 


CPSER(J) [C(.J)-C(J-1)  ] + 


+ SVPSER(J) [C(J)-C(J-1) ] + OCPSER'J)*C(J) 

+ CPSP(.T)  (Y(J)-Y(J-l)  ] + + SVPSP(J)  [Y(J)-Y(J-1)  ]” 

+ CIC(.T)  *Y( J)  + CIPC(.T)  + [URC ( J)  + UTC(J) ]RBAR(J) J , 

where  the  symbols  are  as  defined  in  Tables  T and  III,  and  the 
[a-b]  ([a-b|  ) indicates  the  maximum  (minimum)  of  (a-b,0)  . 


(5) 


Now  the  heuristic  algorithm  present  worth  is  taken  to  be 


PR- WORTH 


+ + C„(.I)  [Y ( J)-Y(  J-l  ) ] + 


, (6) 


where  and  C.n  are  given  by 


C = CPSEK  + OCPSER n.1+.R)-LG-tR.) . . li 1~|  - SVPSER  [" Vt+tI 

1 L ro+r)k-i+1  J L(f+R)K_1+1  J 

(7) 

r(1+R)[(i+R),K"1+1-  ul  _ svps^  r — ] — 1 . 

L Rd+R)™  J L(l+R)K"i+1J 


C,,  = CPSP  + CIC 


First,  the  and  are  computed,  where  is  a function  of 

the  purchase  cost,  operating  cost,  and  salvage  value  of  a repair  channel 
and  is  a function  of  the  purchase  cost,  carrying  cost,  and  salvage 

value  of  a spare  as  given  by  Equation  (7)  for  the  dynamic  model.  The 
first  bracket  term  brings  the  annual  costs,  OCPSER  ind  CIC,  to  a beginning 
of  year  i equivalent  cost,  while  the  second  bracket  term  brings  the  sal- 
vage value  to  a beginning  of  year  i equivalent  term;  that  is,  the  bracket 
terms  are  the  present  worth  factors  for  a beginning  of  year  series  payment 
and  end  of  horizon  payment,  respectively.  Note  that  the  algorithm  assumes 
that  if  a spare  or  repair  channel  is  purchased  in  year  i,  the  annual  costs 
at  year  i values  are  incurred  through  the  end  of  the  horizon,  even  if 
removal  occurs  sooner. 


- 31  - 


The  algorithm  forms  a ratio  (call  A)  of  C^/C^  or  ^/C  , de- 

pending on  the  relative  magnitudes  in  such  a way  that  the  ratio  is  > 1 . 
Then  given  a pair  of  values  C,Y  (to  start  year  i,  | and  Y^_^  are 

used)  the  availability  is  computed.  If  it  is  below  the  desired  level  and 
if,  for  example,  A = C^/C ^ , then  for  an  equal  dollar  expenditure  A re- 
pair channels  or  one  spare  can  be  added.  Availability  is  calculated  for 
both  cases  (adding  A repair  channels  or  one  spare)  and  the  case  yielding 
the  higher  availability  becomes  the  new  C,Y  pair.  The  algorithm  repeats 
until  the  desired  availability  is  met.  Upon  exceeding  the  desired  avail- 
ability, a backoff  procedure  is  utilized.  If  feasibility  was  reached  by 
adding  A channels,  the  algorithm  first  attempts  to  remove  a spare  and 
then  channels  are  removed  one  at  a time  to  see  if  a cheaper  solution 
exists  near  the  boundary.  If  feasibility  was  reached  by  adding  a spare, 
again  one-au-a-time  removal  of  channels  is  tried.  Had  A = » t'ie 

words  channel  and  spare  would  be  reversed  in  describing  the  algorithm. 

When  the  initial  values  of  C and  Y for  year  i exceed  the  avail- 
ability desired,  the  algorithm  immediately  goes  into  a backoff  mode,  try- 
ing to  remove  spares  and  channels  one  at  a time,  starting  with  the  more 
expensive  (larger  value)  first. 

The  algorithm  uses  only  and  C.,  . The  other  costs  (URC,  UTC, 

CIPC)  are  not  used  in  the  algorithm  but  are  considered  in  the  true  cost 
calculations.  The  costs  inside  the  braces  in  Equations  (5)  and  (6)  are 
what  is  given  as  TRUE-COST  and  COST,  respectively,  for  each  year  in  the 
output . 


In  the  static  mode  the  algorithm  works  in  the  same  way,  except 
the  functions  C^  and  C.^  are  changed  to  reflect  all  costs  as  equiva- 
lent uniform  series  end  of  period  expenditures  over  the  system  life. 
Thus,  the  purchase  costs  and  salvage  values  are  multiplied  by  sinking 
fund  and  capital  recovery  factors,  and  the  yearly  operating  costs  asso- 
ciated with  spares  and  channels  which  are  assumed  beginning  of 
period  expenditures  are  multiplied  by  ( 1+R)  . Hence, 
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The  costs  URC,  UTC,  and  CIPC  are  also  assumed  year  beginning  costs 
and  are  multiplied  by  (1+R)  to  bring  them  to  year-end  expenditures,  and 
are  incorporated  into  the  TRUE-COST  calculation  by  adding  (URC  + UTC)  • 
(1+R)R  and  CIPC*(1+R)  to  C^  x (C)  + Cf>  x (Y)  . This  is  then  the  value 


which  shows  as  both  TRUE-COST  and  TRUE-PW  on  the  output,  TRUE-PW  (as  well 
as  PR-WORTH,  which  equals  COST)  being  redundant  ir  the  static  mode. 


A sample  of  intermediate  output  is  shown  in  Figure  7 for  the 
first  year  of  Sample  Run  1.  Shown  are  the  failure  rate  for  year  i (RLAM), 
average  population  failure  rate  for  year  i (AMTBR) , average  turn-around 
(repair)  time  (ST),  availability  for  the  particular  combination  of  C and 
Y,  average  queue  size  at  repair  depot  (EQ),  and  average  number  of  units 
in  repair  (I,). 
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